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This paper shows that under certain symmetry conditions the probability of remaining in the initial 
state (the probabihty of no transition) in a chainwise-connected multistate system driven by two or 
more delayed laser pulses does not depend on the pulse order. 
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The process of stimulated Raman adiabatic passage 
(STIRAP) has received a great deal of attention in the 
past decade ||,|^ because of its potential for efficient and 
robust population transfer between two states ipi and ■03 
via an intermediate state tfj2- STIRAP uses two delayed 
laser pulses, a pump pulse Qp{t) linking states V'l and 
■02 and a Stokes pulse fls{t) linking states ■02 and ■03. By 
applying the Stokes pulse before the pump pulse (coun- 
terintuitive order) and maintaining adiabatic-evolution 
conditions and two-photon resonance between states ■0i 
and ^3 , one ensures complete and smooth transfer of pop- 
ulation from ^01 to ^03, regardless of whether the interme- 
diate state is on or off single-photon resonance. Apply- 
ing the two pulses in the intuitive order [ilp{t) before 
n5(i)] leads to oscillations in the on-resonance case and 
to STIRAP-like transfer in the off-resonance case. The 
success of STIRAP has prompted its extension to mul- 
tistate chainwise-connected systems y-^, where a simi- 
lar distinction between the intuitive and counterintuitive 
pulse orders exists. 

In view of the great difference in the final-state popula- 
tion for the two pulse orders, surprisingly, the initial-state 
population has been found to be the same for both orders 
in the three-state case, provided the Hamiltonian has a 
certain symmetry |^. The present paper extends this 
result to multistate chains. Thus it establishes another 
similarity between three-state and multistate systems. 

The time evolutions of the probability amplitudes 
c{t) — [ci(t),C2(t), . . . ,CAr(t)] of the TV states satisfy 
the Schrodinger equation (in units h= 1) |10|, 



ic{t) =H(i)c(i). 



(1) 



In the rotating-wave approximation the Hamiltonian of 
the multistate chain is given by the tridiagonal matrix 
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The system is supposed to have A^ = 2?t, -|- 1 states and 
the Rabi frequencies r2jj+i(i) obey the relations 



O (f\- i O^pW' J odd, 
^h,j+im-<^^^ns{t), J even, 

np{t) = Qofit~T), 



(3a) 

(3b) 
(3c) 
(3d) 



and f{—x) = f{x). The functions np{t) and fls{t) de- 
scribe the envelopes of the two pulses, 2t is the pulse 
delay, flo is an appropriate unit of Rabi frequency, and 
the (constant) relative coupling strengths ^j are propor- 
tional to the corresponding Clebsch-Gordan coefficients. 
The detunings are supposed to obey the relations 



Aj{t) = Ajv 
A,(t) 



-lit), 



A,(-i), 



(i = 2,3,...,n + l). 



(4a) 
(4b) 



For example, Eqs. (^ and (|4|) apply to 2J + 1-state 
systems (J integer), formed by the sublevels in J ^ J or 
J ^ J — 1 transition, coupled by two laser pulses fip {t) 
and fls{t) with cr+ and a^ polarizations |0-p|. 

I shall show that when conditions (0) and (m are satis- 
fied the probability of remaining in the initial state (the 
probability of no transition) does not depend on the pulse 
order, i.e., it is invariant upon the interchange of flp{t) 
and Qsit)- Since the Qp <=^ fls swap is equivalent to the 
index change j ^ N -{- 1 — j in 'H.{t), the flp <=^ ils in- 
variance of the population of the initial state ipj is equiv- 
alent to the assertion that for a given pulse order, the 
probability of remaining in state ipj, provided the sys- 
tem is initially in state ijjj, is equal to the probability of 
remaining in state ■0jv+i-j, provided the system is ini- 
tially in state ipN+i-j- In terms of the transition matrix 
U(+oo, —00), defined by c(+oo) = U(+oo, — cx))c(— cx)), 
this invariance means that for any j = l,2,...,n + l. 



[/,-j( + 00, -00) = UN-f-l-j,N+l~j{ + 00, -00). 



(5) 



The proof of Eq. (||) is carried out in several steps. The 
first step is to show that the eigenvalues and the eigen- 
states of H(i) have certain symmetric properties. These 
properties lead to symmetries of the Hamiltonian in the 
adiabatic basis, which determine certain symmetries of 
the adiabatic transition matrix, which in turn lead to 
the property (^ of the diabatic transition matrix. 

It follows from Eqs. (||) and (||) that the flp ^ fls 
swap is equivalent to time reversal in H(i), 



i7p(t) ^ ^s{t) is equivalent to 



t 



-t. 



(6) 



Hence, since the J7p <=^ f^s swap does not change the 
eigenvalues of the Hamiltonian, H(— i) has the same 
eigenvalues as H(i). The eigenvalues Xj{t) of H(/:) are 
therefore even functions of time, 



A,(-i) = A,(i), (j = l,2,...,iV). 



(7) 



Since H(i) is real and symmetric, its eigenvalues are real 
and its eigenstates can be chosen real too. The com- 
ponents of the eigenstates (the adiabatic states) w-'(t) — 
[w{{t) , W2{t) , . . . , w]^(i)]'^ are expressed in terms of wi(t) 
(for simplicity, the label j is omitted for the moment) as 



W2it) A(i) 



.92 (i), 



wi{t) Ci^p(i) 

wsit) _ A(t) [A(t) - A2(i)] - ei^Ut) _ 



wiit) 



^i^2npims{t) 



53 (i), 



and in terms of wn (t) as 
WN-l{t) _ Xjt) 



92i-t), 



WN-2{t) _ A(i) [A(i) - A2{t)] - ei^lit) _ 



WN{t) 



ei6^p(t)f^s(i) 



= 93{-t), 



Generally, one can write Wk(t)/wi{t) = gk{t) and 
WN+i-kit)/wN(t) = gk{—t). For fc = n + 1, one finds 
gn+i{—t)w]S!{t) = gn+i{t)wi{t). It follows that 

WN+i-k{t) gki-t) gn+i{t) 



Wkit) gk{t) g„+i(-i)' 

for any k = 1, 2, . . . , n + 1. Hence 

wiit) ^ gn+ii-t)/u{t), 
W2{t) = g2{t)gn+i{-t)/v{t), 

' ' ' 1 

Wn+iit) = g„+i{t)gn+i{-t)/u(t), 

' ' ' ; 

■WN-iit) = gn+i{t)g2{-t)/i^{t), 

WNit) = gn+l{t)/v{t)- 



(8) 



The normalization factor v(t) is obviously invariant upon 
time reversal, which means that v{~t) — i'{t). Equations 
(^, which are valid for g^+i(i) f^ (case I), lead to the 
relation (with the label j restored) 



i.(-t) = w]^+i_fe(i), (easel). 



(9a) 



with fe = l,2,...,n-|-l. 

If gJIYi(^) ~ ^ (case II) for a certain \m{t), we have 
w™^^(i) = and w"+2(-^) = -w'"(^), which leads to 



<(-i) = -^«5^Vi-fe(t), (case II), 



(9b) 



with fc = l,2,...,n+l. Such a case arises for the zero- 
eigenvalue eigenstate in systems with N = 3,7,11,... 
states and zero detunings. 

The symmetry relations (^ for the adiabatic states 
determine certain symmetries of the Hamiltonian in the 
adiabatic basis. The transformation from the original 
(diabatic) basis to the adiabatic basis, c(i) = W(i)a(i), 
is carried out by the orthogonal matrix W(t), whose 
columns are the normalized eigenvectors w-'(i). Here 
a(t) = [ai{t),a2{t), . ■ . ,aN{t)] is the column-vector of 
the adiabatic probability amplitudes. The Schrodinger 
equation in the adiabatic basis reads 



where H"(i) : 

radb 
Trnonadb 



ik{t) =H''(i)a(i), 
H^'^'^(i) + H"°"^'i'^(i) with 



H^a'^(i) =W^(t)H(i)W(i), 
(i) = -iW^(t)W(t). 



(10) 



(11a) 
(lib) 



The adiabatic part H'''^'^(i) is a diagonal matrix contain- 
ing the eigenvalues Aj [t) of H(i) on the main diagonal. 
The nonadiabatic part H"°"^'^'^(i) has zeros on the main 
diagonal, while the off-diagonal elements are equal to the 
nonadiabatic couphngs H^'^'^^{t) = -iw^{t) ■ w''(t). It 



is readily seen from Eq. (fed) that the nonadiabatic cou- 
pling between two case-I adiabatic states w-' (i) and w'^ (t) 
is an odd function of time. Really, 



N 



H^on.6^{^t)^^Y,^i{~tH{-t) 



1=1 

N 



= iJ2^N+l-li*)^N+l-lit) 



1=1 



= -iJj^r (0, (easel- easel). (12a) 

The nonadiabatic coupling between a ease-I eigenstate 
w^{t) and a case-II eigenstate w'"(t) is an even function. 



H^Z {-t) - ^"r (0, (case I • case II). (12b) 

The symmetry of H''(t) determines a certain sym- 
metry of the adiabatic transition matrix U°(-|-oo, — oo), 
defined as a(+oo) = U°(-|-cw, — cx3)a(— cx)). In order 
to find it, I introduce the evolution matrix G(t, 0) via 
a(t) = G(i, O)a(O). Evidently, the first column of G(i, 0) 
is the solution of Eq. (M) for the initial condition a(0) — 

rp ^ 

(1,0,0,. ..,0) , the second column is the solution for 
a(0) = (0, 1, 0, . . . ,0) , and so on. When all nqnadi 



abatic couplings are odd functions of time [Eq. (12a)], 
time reversal in Eq. ( |l0|) is equivalent to complex conju- 
gation of a(i) (case A) . When a case-II eigenstate vif™ [t) 
exists [then the nonadiabatic couplings involving it are 
even functions, Eq. ( |12b| )], time reversal in Eq. (Ilfl) is 
equivalent to complex conjugation of a(i) and change of 
sign of am{t) (case B). This means that 



G(-t,0) 



G*(t,0), (case A), 
IG*(t,0)I, (caseB), 



(13) 



where I is a diagonal matrix with units on its diagonal, 
except the {m, m)-th element which is —1. It follows from 
Eq. (O) and the unitarity of G that 



U"(+oo,-oo) 



G(+oo,0)G(0,-oo) 

G(+cx),0)G^(-oo,0) 

G(+cx),0)G'^(+cx),0), (case A), 
G(+(X),0)IG^(+cx),0)I, (caseB). 



Hence 

[U"(+ 



"^)Y 



U''(+oo, — oo), (case A), 
IU"(+oo,-oo)I, (caseB). 



(14) 



The transition matrices in the diabatic and adiabatic 
bases are related by 

U(+oo, -oo) = W(+oo)U"(+oo, -cx))W'^(-cx)). 

Then one finds from Eqs. (g) and ( [l4[) that in both cases 
A and B, 



UN+l~j,N+l-ji + 00, -oo) 

N 

N 

= Y. '^j^(-oo)C/ffc(+oo,-oo)^(+oo) 

k,l = l 

= Ujj{+oo,~oo). 

This completes the proof. 

It should be emphasized that the pulse-order invari- 
ance applies to the population of the initial state only, 
while the populations of all other states depend on the 
pulse order. This is clearly demonstrated in Fig. m where 
the initial-state population Pi and the final-state popula- 
tion P5 are plotted against the pulse delay r in the case 
of a five-state system, initially in state ipi. The figure 
shows that P5 behaves very similarly to STIRAP with a 
broad plateau of high transfer efficiency for r > and os- 
cillations for T < [QplO]. In contrast. Pi is a symmetric 
function of r, as follows from the above results. 

Finally, the pulse-order invariance of the initial-state 
population has been derived without the assumption of 
adiabatic evolution. Hence it applies to the general nona- 
diabatic case, as long as the pulse duration is long enough 
to validate the rotating-wave approximation. 

This work has been supported financially by the 
Academy of Finland. 
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FIG. 1. The initial-state population Pi and the fi- 

nal-state population P5 for a five-state system, initially in 
state tpi, plotted against the pulse delay r in the reso- 
nance case (A2 — A3 = A4 — 0). The Rabi fre- 
quencies of the two pulses are given by Eqs. dd) and 
have Gaussian shapes, Qp{t) = r2oexp[— (f — t)'^/T^] and 
ns{t) = noexp[-(i -H r)Vr^], with ^1=^4 = \/T/3, 



6=^3 



''1/2, and QqT = 30. 
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